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NOTATION 

Are* 

Axisymmetric  dispiacefflent  thickness,  Equation  [12] 

Constant  in  Mangier**  transformations,  Equation  [42] 

Drag  coefficient,  D/—p  U2 
2 

Total  drag 
Sktn*fricUon  drag 

Pressure  drag 

Subscript  for  quantities  at  tail  of  body 

Subscript  for  quantities  at  point  where  S <<  no  longer  holds 
Two-dimensional  shape  parameter,  S*/6 
Axisymmetric  shape  parameter,  A*/D 
Exponent  in  Equation  [76],  (1  -f  m)  (ff  2) 

Length  of  body 

Axial  distance  along  body  of  revolution 

Exponent  of  power  law  for  local  coefficient  of  skin  friction. 
Equation  [74] 

Subscript  for  quantities  at  neutral  stability  point 

Exponent  of  power  law  for  velocity  pronies.  Equation  [85] 

Pressure  in  boundary  layer 

Pressure  of  undisturbed  incoming  stream 

Exponent  in  Equation  [92] 

Velocity  of  fluid  in  ^-direction 

Turbulent  fluctuation  in 

Body  Reynolds  number,  LVjv 

Boundary-layer  Reynolds  number,  OV/v 

Transverse  radius  from  axis  of  body 

Transverse  radius  to  surface  of  body  of  revolution 

Radius  to  equivalent  body  with  added  displacement  thickness 

Subscript  for  quantities  at  transition 

Subscript  for  quantities  at  transition  point  induced  by  Iree- 
stream  turbulence 
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Subscript  for  (pisatitlM  at  sslf*exoit«d  transition  point 

Vnlooity  at  outar  adgo  of  boundary  layer 

Velocity  of  undisturbed  incomiiig  stream 

Velocity  of  fluid  parallel  to  surface  of  body  in  »>direction 
(mean  velocity  for  turbulent  flows) 

Turbulent  fluctuation  in  w 

Velocity  of  fluid  normal  to  surface  of  body  (mean  velocity 
for  turbulent  flows) 

Turbulent  fluctuation  in  v 

Arc  length  along  meridian  profile 

Distance  normal  to  body  surface 


Angle  between  tangent  to  surface  and  axis  of  body 
Thickness  of  boundary  layer 
Two*dimenaional  displacement  thickness 


Constant  in  power  law  for  local  coefficient  of  skin  friction 
in  Equation  [74) 


Two'dimensional  momentum  thickness 


Displacement  area, 


Viscosity  of  fluid 
Kinematic  viscosity  of  fluid 
Density  of  fluid 
Normal-stress  term, 


Shearing  stress,  g- 


Shearing  stress  at  body  surface 

^ Angle  between  a meridian  plane  and  the  reference  meridian  plane 

n Momentum  area. 


0 


Momentum  area  of  wake  at  infinity  downstream 

Subscript  for  quantities  at  xero  pressure  gradient  on  flat  plates 
Tilde  over  equivalent  two-dimensional  quantities 


ABSTRACT 


A procedure  is  presented  for  oelculnting  the  viscous  dreg  oC  bodies  of 
revolution  in  sxisl  motion  from  boundsry*lsyer  theory.  Rspid  spproximste 
methods  are  developed  for  computing  the  growth  of  the  Isminsr  end  turbulent 
boundary  layers.  A new  empirical  criterion  is  given  for  locating  the  position  of 
self«excited  transition  associated  with  low*turbulence  flows. 

INTRODUCTION 

A completely  immersed  body  moving  rectilinearly  with  uniform  velocity  in  an  infinite 
fluid  at  rest  experiences  a resisting  force  which  may  be  termed  viscous  drag  as  it  results  pri* 
marily  from  the  viscous  properties  of  the  fluid.  In  designing  low*drag  bodies  there  often  a> 
rises  the  need  for  calculating  the  viscous  drag  of  streamlined  bodies  of  revolution  in  axial 
motion  when  considering  various  proposed  shapes.  For  bodies  moving  at  high  Reynolds  num- 
bers, which  are  of  great  technical  importance,  the  viscous  drsg  of  streamlined  bodies  of  revo- 
lution is  readily  amenable  to  analytical  treatment  on  the  basis  of  the  boundary-layer  concept. 

Historically,  the  theoretical  analysis  of  the  drag  of  bodies  in  uniform  motion  by  as- 
suming an  ideal  (non-viscous)  fluid  gave  the  fruitless  result  of  teto  drag  for  all  bodies,  the 
classical  D'Alembert  paradox.  At  the  other  extreme,  the  theoretical  analysis  of  drag  by  apply- 
ing the  complete  set  of  Navier-Stokes  equations  of  motion  for  the  flow  of  a viscous  fluid  led, 
in  general,  to  mathematical  difficulties  which  were  virtually  unresolvable  owing  to  the  com- 
plicated non-linear  nature  of  these  equations.  For  bodies  moving  at  high  Reynolds  numbers, 
however,  the  flow  is  virtually  that  of  an  ideal  fluid  except  in  a thin  boundary  layer  next  to  the 
body  where  substantial  viscous  forces  are  produced  by  the  rapid  drop  in  velocity  to  aero  at  the 
boc^  surface.  Accordingly,  by  omsidering  the  viscous  flow  confined  to  the  boundary  layer, 
Prandtl  was  able  to  derive  the  simpler  boundary-layer  equations  of  motion  from  the  Navier- 
Stokes  equations. 

The  principal  purpose  of  this  report  is  to  describe  methods  of  solving  the  boun&ry- 
layer  equations  of  motion  to  arrive  at  the  viscous  drag  of  bodies  of  revrdution  of  arbitrary 
shape  in  uniform  axial  motion.  The  study  is  restricted  to  hydraulically  or  aerodynamically 
smooth  streamlined  bodies  in  incompressible  flow.  A streamlined  body  may  be  defined  as  one 
without  appreciable  separation  of  flow  from  its  surface  and  consequently  with  small  pressure 
drag  resulting  from  the  generation  of  separation  eddies.  It  is  to  be  noted,  however,  that  some 
pressure  drag  is  still  present  in  the  viscous  drag  of  even  perfectly  streamlined  shapes  owing 
to  the  effect  of  the  boundary  layer  in  displacing  the  main  flow  outward,  especially  near  the 
tail. 

The  calculation  of  the  viscous  ckag  of  a body  of  revolution  requires  a detailed  analysis 
of  the  development  of  each  phase  of  the  boundary  layer  from  its  origin  on  the  nose  of  the  body 
to  its  final  phase  as  the  frictional  wake  far  downstream.  In  the  downstream  direction  the 


boundary  layer  may  consist  successively  oft  a laminar  boundary  layer,  a transition  tone  from 
laminar  to  turbulent  flow,  a turbulent  boundary  layer,and  a McUonal  wake. 

In  the  laminar  boundary  layer  an  approximate  method  involving  a simple  quadrature  is 
derived  for  the  rapid  calculation  of  the  changes  in  momentum.  The  derivation  consists  of  an 
extension  to  axisymmetric  flows  past  bodies  of  revolution  of  a method  of  successive  approx* 
imation  introduced  by  Shvets*  for  two-dimensional  laminar  boundary  layers. 

New  empirical  criteria  are  presented  for  locating  the  position  of  transition  for  either 
low-turbulence  or  turbulent  free-streams  from  the  position  of  neutral  stability.  The  position 
of  so-called  Self-excited  transition  occurring  in  low-turbulence  streams  or  under  flight  con- 
ditions is  based  on  the  average  pressure  gradient  from  the  position  of  neutral  stability  to  that 
of  transition.  Although  the  test  data  are  (or  two-dimensional  flows,  the  criterion  is  extended 
to  axisymmetric  flows  past  bodies  of  revolution  by  means  of  Mangier*  s transformations.  An 
approximate  criterion  for  estimating  the  position  of  transition  on  a body  in  a turbulent  free- 
stream  is  established  on  the  basis  of  measured  positions  of  transition  for  flat  plates  in  free 
streams  with  various  degrees  of  turbulence. 

The  analysis  of  the  axisymmetric  turbulent  boundary  layer  on  a body  of  revolution  is 
divided  into  that  for  the  main  portion  of  the  body,  where  the  boundary  layer  is  relatively  thin 
compared  to  the  radius  of  the  body,  and  that  for  the  tail  portion,  where  the  boundary  layer  is 
relatively  thick.  The  momentum  changes  in  the  thin  boundary  layer  may  be  calculated  by  a 
rapid  method  involving  a simple  quadrature  wherein  a power-law  relation  for  the  skin  friction 
of  flat  plates  is  incorporated.  Flat-plate  values  for  skin  friction  are  deemed  reliable  even 
where  the  local  skin  friction  is  diminishing  in  an  adverse  pressure  gradient  owing  to  the  com- 
pensating effect  of  the  Reynolds  normal-stress  term.  The  turbulent  boundary  layer  on  the  tail 
is  analyzed  by  means  of  appropriate  linear  simplifications  which  give  a rapid  method  for  cal- 
culating the  momentum  changes.  An  expression  is  derived  for  the  change  in  momentum  pro- 
duced by  the  pressure  difference  in  the  wake  at  the  tail  and  in  the  wake  far  downstream  in  ac- 
cordance with  a method  presented  by  Young^  wherein,  however,  a nrore  general  relationship  is 
employed  for  the  variation  of  the  shape  parameter  of  the  velocity  pronie. 

For  quick  reference  the  various  steps  involved  in  calculating  the  development  of  the 
boundary  layer  on  a body  of  revolution  are  summarized  at  the  end  of  this  report. 

GENERAL  CONSIDERATIONS 
AXISYMMETRIC  BOUNDARY-LAYER  FLOW 

The  main  elements  comprising  axisymmetric  flow  past  a body  of  revolution  at  high 
Reynolds  numbers  are  shown  in  Figure  1 for  the  meridian  plane.  Two  principal  regions  of 
flow  are  indicated:  the  boundary  layer  next  to  the  body  with  viscous  flow  and  the  region 


*Rafar«ncaa  are  listed  on  page  29. 
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extwnal  to  the  boundary  layer  with  essentially  potential  (non*viscous)  How.  The  pressure  dis 
tribution  along  the  body  is  mainly  determined  by  the  potential  Row  with  only  a slight  modiO* 
cation  arising  from  the  boundary-layer  flow  except  right  at  the  tail  where  the  after  stagnation 
pressure  of  potential  flows  is  absent.  The  pressure  distribution  on  the  body  is  important  in 
determining  the  boundary* layer  flow  inasmuch  as  it  is  one  of  the  most  important  factors  gov- 
erning the  growth  and  development  of  the  boundary  layer.  A negative  pressure  gradient  is 
termed  favorable  and  a positive  pressure  gradient  adverse  in  connection  with  preventing  sepa- 
ration of  the  boundary-layer  flow  from  the  body  surface. 

Starting  as  laminar  flow  at  the  stagnation  point  on  the  nose,  the  boundary-layer  flow  de- 
velops instability  and  undergoes  transition  to  turbulent  flow  at  some  position  downstream  on 
the  body.  The  boundary-layer  flow  usually  continues  turbulent  for  the  remaining  after  portion 
of  the  body  and  leaves  the  tail  as  the  frictional  wake  which  extends  indefinitely  downstream. 

VISCOUS  DRAG 

The  viscous  drag  of  a bod^  is  generally  derivable  from  the  boundary-layer  flow  either 
on  the  basis  of  the  local  forces  acting  on  the  surface  of  the  body  or  on  the  basis  of  the  ve- 
locity profile  of  the  wake  far  downstream.  The  local  hydrodynamic  force  on  a unit  of  surface 
area  is  resolvable  into  a surface  shearing  stress  or  local  skin  friction  tangent  to  the  body 
surfoce  and  a pressure  p normal  to  the  surface.  The  summation  over  the  whole  body  surface 
of  the  axial  components  of  the  local  skin  friction  and  of  the  pressure  gives,  respectively,  the 
skin-friction  drag  Dj  and  the  pressure  drag  which  for  a body  of  revolution  in  axisym'metric 
flow  become 
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Figure  1 - Typical  Boundary  Layer  Around  a Bo<ty  of  Revolution  in  a Meridian  Plane 
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where  Is  the  radius  from  the  axis  to  the  body  surface, 

a is  the  angle  between  the  tangent  to  the  meridian  profile  and  the  axis  of  the  body, 
e is  the  arc  length  along  the  meridian' profile,  and 
e,  is  the  total  arc  length  of  the  body  from  nose  to  tail. 

The  sum  of  the  two  drage  then  constitutes  the  total  viscous  drag  D or 

^ [31 

The  alternate  procedure  involving  the  velocity  profile  of  the  wake,  which  is  better 
suited  to  the  methods  of  this  report,  is  to  determine  the  total  drag,  without  need  of  resolution 
into  skin«friction  and  pressure  drag,  by  considering  the  net  rate  of  loss  of  momentum  of  the 
flow  of  the  entire  stream.  The  analytical  procedure  is  to  apply  the  momentum  theorem  of 
hydrodynamics  to  a control  surfrice  enclosing  a region  about  the  body  with  dimensions  suf* 
ficiently  large  to  have  substantially  undisturbed  free« stream  pressure  at  its  periphery.  The 
total  drag  of  the  body  which  is  the  net  rate  of  loss  of  momentum  in  the  axial  direction  is  then 
given  by 


D»2mpj  u{V^-u)rdr 


where  is  the  velocity  of  the  incoming  undisturbed  stream  ahead  of  the  body, 
u is  the  velocity  in  the  wake  far  downstream, 
r is  the  radial  distance  from  the  axis,  and 
p is  the  mass  density  of  the  fluid. 

In  terms  of  the  momentum  area  of  the  wake  far  downsUeam 


the  drag  is 


"-•1 


D-2np  uj 


The  drag  coefficient  Cq  based  on  some  appropriate  reference  area  A is 

Cq ^ m 

Inasmuch  as  the  momentum  area  of  the  wake  is  the  final  stage  of  the  development 
of  the  boundarylayer  flow  from  its  incepUon  on  the  nose  of  the  body,  it  is  necessary  to  cal* 
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cuUte  the  momentum  area  of  the  boundary*layer  flow  for  each  of  its  stages  in  order  to  detei^ 
mine  the  resulting  drag  of  the  body. 


PRESSURE  DISTRIBUTION  ALONG  THE  BODY 


Before  any  boundary*layer  calculations  can  be  performed  for  the  body  of  revolution 
whose  drag  is  to  be  determined,  it  is  necessary  to  have  on  hand  the  distribution  of  pressure  p 
along  the  body,  which  in  non>dimensional  terms  is  usually  presented  as 


P - A. 
2 “ 


[81 


Here  p^  is  the  pressure  of  the  undisturbed  stream  far  ahead,  and  U is  the  velocity  at  the  outer 
edge  of  the  boundary  layer. 

Where  no  experimental  values  of  pressure  distribution  are  available,  recourse  can  be 
had  to  methods  for  calculating  the  potential  flow  past  the  body  which  involve  the  solution  of 
the  Laplace  equation  for  arbitrary  boundary  conditions.  Such  methods,  which  are  reviewed 
briefly  and  evaluated  in  Reference  3,  are,  in  general,  numerically  arduous  and  difficult  to 
apply  to  bodies  of  arbitrary  shape.  Recently,  Landweber^  developed  an  accurate  method,  well 
suited  to  automatic  calculating  machines,  in  which  a special  iteration  formula  is  employed  in 
solving  the  resulting  Fredholm  integral  equation  of  the  first  kind.  A faster  method,  giving, 
however,  only  approximate  results,  is  that  of  Young  and  Owen^  which  involves  interpolation 
among  tabulated  values  of  Legendre  polynomials. 

For  the  original  profile  of  the  body  of  revolution  the  calculated  pressure  distribution 
for  potential  flow  agrees  closely  with  measured  values  over  most  of  the  length  of  the  body, 
the  greatest  discrepancy  appearing  neat' the  tail  because  of  the  displacement  effect  of  the 
boundary  layer.  A closer  result  can  be  obtained,  however,  by  repeating  the  potential*flow  cal- 
culation for  a somewhat  altered  body  consisting  of  the  original  contour  and  an  added  thickness 
based  on  the  displacement  effect  of  the  boundary  layer.  The  added  displacement  thickness  in 
the  y*direction  normal  to  the  surface  for  a body  of  revolution  in  axisymmetric  flow,  as  shown 
in  Figure  2,  is  obtained  by  equating  the  total  flow  retarded  in  the  boundary  layer  of  thickness 
5 to  the  amount  subtracted  from  the  potential  flow  of  thickness  a*  or 


or 


o* 

J 2ff(f/-0)rdy 
0 


J 2n{V  -u)  rdy 


a* 


8 


\rdy.  J - A* 


[9] 
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Figure  2 • Equivalent  Body  of  Revolution  with  Displacement  Thickness 
where  A*  is  the  displacement  area.  Substituting 


till 

tl2] 

1131 


which  is  tlie  displacement  thickness  for  two-dimensional  flow.  In  the  wake  where  - 0, 
Equation  [12]  reduces  to 


[141 
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It  is  convenient  to  define  n rsdins  r*  to  the  surface  of  the  equivalent  body,  shown  in 
Figure  2, 

f*  - + a cos  a [15] 

Then  from  [12] 

f*  - + 2 A*  coo  a [16] 


AXISYMMETRIC  LAMINAR  BOUNDARY  LAYER 

The  solution  of  laminar  boundary  layers  in  either  two*dimensionai  or  axisymmetiic 
flows  has  been  the  subject  of  numerous  mathematical  investigations  owing  to  the  interesting 
characteristics  of  the  resulting  equations  of  motion.  The  considerable  literature  that  has 
developed  contains  methods  of  various  degrees  of  complexity  and  precision  which  are  given  in 
the  summaries  by  Goldstein^  and  more  recently  in  the  AVA  Monographs^  and  by  SchUchting.^ 
For  many  drag  calculations  very  precise  solutions  of  the  boundary  layer,  requiring  extensive 
numerical  work,  are  usually  not  warranted,  especially  where  the  laminar  boundary  layer  con- 
stitutes a small  part  of  the  whole  boundary  layer.  A simple  approximate  formula,  well  adapted 
to  drag  calculations,  will  be  devised  by  an  extension  to  axisymmetric  flow  of  a method  of  suc- 
cessive approximation  introduced  by  Shvets^  for  two-dimensional  flows. 

Owing  to  the  thinness  of  the  laminar  boundary  layer  on  the  forward  part  of  the  body, 

5 « r^ , the  equations  of  motion  for  steady  axisymmetric  boundary-layer  flow  past  a body  of 
revolution  with  negligible  longitudinal  curvature  (Reference  5)  reduce  to 


--i  il  + v ^ 

dy  p d9 


o-i£ 

9y 


and  the  equation  of  contin  a^y  reduces  to 


-f  JL 
V ^ 


[17] 


[18] 


Here  u and  v are  the  «-  and  y-components  of  the  boundary-layer  velocity  respectively  parallel 
and  normal  to  the  surface  of  the  body  in  the  meridian  plane,  p is  the  pressure  in  the  boundary 
layer,  v is  the  kinematic  viscosity  of  the  fluid,  and  S is  the  thickness  of  the  boundary  layer 
in  the  y-direction.  These  equations  have  been  shown*  to  remain  applicable  at  the  forward 
stagnation  point  for  bodies  with  blunt  noses  even  thdugh  both  0.  The  following 

boundary  conditions  are  to  be  satisfied  by  the  boundary-layer  equations: 


8 
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u - 0 at  y - 0 
u»V  at  y mi 


tl9l 


The  pressure  p and  the  velocity  U at  the  outer  edge  of  the  boundary  layer  are  related  by 
Bemouilli's  equation  for  potential  flow 

f - constant  [20] 


or 


-i$L-  u HL 

p Mt  dm 


[21] 


Combining  the  equation  of  motion,  [17],  the  continuity  equation,  [18],  and  the  differential 
Bemouilli  equation,  [21],  produces 


dy2 


[22] 


A first  approximation  is  to  let 


[23] 


Integrating  [23]  twice  and  utilizing  boundary  conditions,  [19],  results  in  a linear  velocity 
ptoftle 


ju 

V 


3L 

8 


[24] 


A second  approximation  is  obtained  by  substituting  the  linear  velocity  profile,  [24], 
into  differential  equation  [22].  Integrating  twice  and  utilizing  boundary  conditions,  [19],  as 
before  gives 


£.iIL 

24  dm 


♦ 


[25] 


Applying  the  additional  boundary  condition  that 


4^.0  at 
^9 


y 


to  the  differentiated  form  of  Equation  [25]  yields 


8 


[26] 


A4£(r 

V dm  ' 


V 


[27] 
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For  A lioMT  v«looity  profiU,  [24], 


4-1 


wh«i«  9 is  the  two>dimen8ionsl  mosMstuin  thiokness 


Substituting  6 in  [27]  gives 


**r 

<(>.»)*,  t iUi. 

dm  C/  4b  ' •*  ' 9 t/ 


Inte^mting  the  linssr  differential  equation,  [SO],  produces 


with  the  stagnation  point  s - 0 as  the  initial  point  of  integration.  An  interesting  feature  of 
this  result  is  its  being  identical  in  both  fom  and  nunerical  constants  to  the  seni-empirical 
result  of  Thwaites*  when  his  twokUmensional  solution  is  extended  to  the  axisynmetric 
case.^^  Thwaites  averaged  known  solutions  of  two^nensiooal  laminsr  boundsry  layers  for 
a variety  of  pressure  gradients. 

Now  the  nx>mentum  area  Q for  axisynunetric  boundary>layor  flow  in  general  may  be  de> 
fined  as 


Sinoe 


r - r 4 y cos  o 


then  in  general 


0 ^ oosa-1 


9 < n < (f,_  + 9)  9 


Hence  for  5 « which  is  the  case  under  consideration 


n - p^9 


[35] 
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A noB*dimeo8ioBal  fomt  of  tho  momwituni  equBtion,  [SI],  oonvonient  for  OBlottlaUon 
putpoaes  ia 


where  I ia  the  axial  diatance  horn  the  aoae, 

dl  cotta dx 

£,  ia  the  length  of  the  body  and  ia  the  body  Reyaolda  number  given  by 

[38] 

Tranaition  to  turbulent  flow  ia  aaaumed  to  occur  inatantaneoualy  at  a tranaition  point 
(1/L),  which  then  becomea  the  upper  limit  of  integration  in  Equation  [86]  when  the  integrataon 
ia  performed  for  the  complete  laminar  boundary  layer. 

TRANSITION 

Aa  the  boundary  layer  thickena  dowiiatream  on  the  body,  the  laminar  flow  tenda  to  be- 
come unatable  and  undergo  tranaition  to  turbulent  flow  under  the  atimulne  of  diaturbanoea  in 
the  flow.  The  tranaition  lone  may  be  conaidered  to  extend  from  the  point  where  the  charao- 
teriatic  ahape  of  the  mean-velocity  profile  of  the  laminar  boundary  layer  begina  to  change  to 
the  point  where  the  characteriatic  ahape  of  the  mean-velocity  profile  of  the  turbulent  boundary 
layer  firat  appeara.  For  moat  drag  calculationa  the  cone  of  tranaition  ia  abort  enough  to  be 
adequately  repreaented  by  a tranaition  point.  The  poeition  of  tranaition  dependa  largely  upon 
the  interaction  of  the  boundary-layer  flow  with  random  diaturbancea  in  the  flow.  Sigpificant 
parametera  of  the  boundary-layer  flow  affecting  the  poaition  of  tranaition  are  the  boundary- 
layer  Reynolda  number  repreaenting  the  ratio  of  inertial  forcea  to  viacoua  forcea,  the  prea- 
aure  gradient  in  the  downatream  direction,  and  the  curvature  of  the  aurface.  The  aource  of 
random  diaturbancea  may  be  the  turbulence  in  the  free  atream,  the  roughneaa  of  the  aurface 
or  noiae  being  tranamitted  through  the  fluid. 

Aa  ahown  theoretically  by  Tollmien  ^ and  Schlichting^  among  othera  (aee  aummary  in 
Reference  12)  and  verified  experimentally  by  Schubauer  and  Skramatad^  and  Liepmann,^^  the 
laminar  boundary  layer  exhibita  atability  characteriatica  which  are  governed  largely  by  the 
boundary-layer  Reynolda  number  and  by  the  preaaure  gradient.  Random  diaturbancea  of  van- 
iahingly  amall  amplitude  have  certain  frequenciea  amplified  and  other  frequenciea  damped  by 
the  laminar  flow  in  Uie  boundary  layer.  The  amplified  fluctuationa  combine  into  regular  wavea 
termed  Tollmien-SchUohting  wavea  which  increaae  in  amplitude  downatream  at  a rate  deter- 
mined by  the  Reynolda  number  and  the  preaaure  gradient  of  the  boundary  layer.  Intermittent 


u 


bursts  of  high'ftoqusBoy  fluotusUons  which  srs  usooistsd  with  intermittent  separntions  of  the 
Inminnr  flow  herald  the  arrival  of  the  more  random  fluotnationa  oharaoteristic  of  fully  devel- 
oped turbulent  flow.  The  transition  Just  deaoribed,  which  depends  on  the  amplifleation  of  vu- 
ishingly  small  disturbances,  may  bs  termed  self-excited  traasitica. 

Disturbances  of  greater  amplitude  arising  from  ftee-stream  turbulence  or  rough  surfaces 
tend  to  hasten  the  decomposition  of  the  laminar  boundary  layer  into  a turbulent  boundary  layer. 
Taylor^’  has  developed  the  concept  of  momentary  separations  arising  from  momentary  adverse 
pressure  gradients  as  the  mechanism  iastigatiag  turbulent  motion.  In  accordance  with  Taylor’s 
analysis  both  the  scale  and  intensity  of  the  free-stream  turbulence  have  been  Aown  experi- 
mentally * to  have  a marked  bearing  on  the  poaitioo  of  transition. 

The  studies  of  Liepmann^^*^^  on  the  stability  of  laminar  boundary  layers  on  curved 
surfaces  have  shown  the  Tollmien-Schlichting  type  of  stability  to  exist  on  surfaces  convex  to 
the  flow  and  the  Gortler  type  of  stability  involving  vortices  to  exist  on  concave  surfaces.  Con- 
vex surfaces  have  a greater  stabilising  effect  and  concave  surfaces  a lesser  stabilising  effect 
than  flat  surfaces. 

The  presence  both  of  large  boundary-layer  Reynolds  numbers  and  of  adverse  pressure 
gradients  tends  to  promote  transition  by  increasing  the  instability  of  the  laminar  boundary 
layer  and  by  accelerating  the  amplification  of  the  ToUmien-Schlichting  waves. 

Under  special  circumstances  transition  is  hastened  when  the  laminar  boundary  layer 
separates  from  the  body  and  reattaches  itself  as  a turbulent  boundary  layer.  Such  separation 
may  be  caused  by  sharp  adverse  pressure  gradients**on  bodies  at  large  angles  of  attack  or  by 
sharp  adverse  pressure  gradients  induced  by  large  single  roughnesses  obstructing  the  flow\^ 

Quantitative  criteria  for  establishing  the  transition  points  on  smooth  bodies  of  revolu- 
tion will  be  considered  here  for  two  technically  impertant  flow  situations  which  are  character- 
ised by  the  absence  or  presence  of  significant  amounts  of  turbulence  in  the  main  flow  stream. 
Zero  or  low-turbulence  condition  exists  in  specislly  constructed  low -turbulence  wind  tunnels 
while  a condition  of  various  degrees  of  turbulence  is  present  in  most  wind  tunnels  and  flow 
facilities.  The  flight  of  aircraft  is  considered  a ease  of  low  turbulence  on  the  basis  of  tests 
by  Jones^^  who  concluded  that  the  scale  of  turbulence  in  the  atmosphere  is  such  as  to  have 
no  effect  on  transition.  A c-milsr  low-turbulence  condition  may  be  assumed  in  the  case  of 
bodies  moving  in  the  depths  of  the  oceans  where  currents  are  absent. 

Transititm  in  the  low-turbulence  case  may  be  considered  to  be  of  the  ToUmien- 
Schlichting  type  wherein  vanishingly  small  disturbances  are  amplified  in  the  boundary  layer 
to  the  castastrophic  point  of  resulting  turbulence  within  the  boundary  layer.  The  stability 
analysis  of  laminar  boundary  layers  shows  the  existence  of  a point  of  neutral  stability  where- 
in the  immediate  neighborhood  upstream  disturbances  of  all  frequencies  are  damped  out. 
Mangier^*  has  prepared  a chart,  partly  reproduced  in  Figure  8,  which  specifies  the  neutral 
stability  point  in  terms  of  the  critical  value  jy  of  boundary-layer  Reynolds  number  Rq 
where 
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— » Prmurt  Crooitm  Poromattr 

Figure  8 - Neutral  Stability  Point  as  FuactioD  of  Preaaure  Qradient  (Mangier,  Ref.  21) 


^ 1391 

as  a function  of  a presaure*flyadient  parameter  ^ ^ . The  curve  ahown,  which  is  baaed  on 

V dm 

the  stability  characteriatica  of  Pohlhausen4ype  velocity  profiles  as  calculated  by  Schlichting 
and  Ulrich^ ^ is  applicable  to  both  two«diinensional  and  axisyimnetric  boundary  layers. 

Unfortunately  no  theoretical  analysis  has  been  developed  for  locating  the  point  of  self* 
excited  transition  from  the  amplification  of  the  Tollmien-Schlichting  waves  past  the  point  of 
neutral  stability.  Experimental  evidence*^  indicates  a pronounced  effect  by  favorable  pressure 
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gradients  in  damping  and  adverse  pressure  gradients  in  amplifying  the  Toiimien*Scblichting 
waves. 

A rough  empirical  criterion  for  locating  the  point  of  self*excited  transitions  is  ascribed 
to  Walz^^  who  suggests  a critical  value  of  the  Reynolds  number  of  transition  as  thrice 
that  for  neutral  stability  or 


[401 

There  follows,  however,  a new  empirical  criterion  which  is  based  on  more  rational 
grounds  and  which  attempts  to  incorporate  the  effect  of  pressure  gradients  in  determining  the 
position  of  self*excitod  transition.  Since  the  transition  point  depends  on  the  cumulative  ef- 
fect of  the  pressure  gradients  from  neutral  stability  to  transition,  it  would  seem  appropriate 
to  use  as  a first  approximation  the  average  pressure  gradient  for  one  of  the  significant  parame- 
ters. Accordingly  two-dimensional  data  on  transition  points  from  tests  on  wing  sections  in 
flight^*^^  and  in  low-turbulence  wind  tunnels^^*^^*^^  were  analyzed  on  this  basis.  As  shown 
in  Figure  4 the  difference  in  Reynolds  numbers  from  the  point  of  neutral  stability  to  transition 

^0  tt  N plotted  against  the  average  pressure  gradient  parameter  ^ over  the 
* * V dx" 

same  region  where 


V dS 


[41J 


The  tildes  over  the  symbols  refer  to  two-dimensional  flows.  Examination  of  the  plotted  data 
in  Figure  4 indicates  a reasonably  consistent  variation  between  the  two  parameters  involved. 

The  conversion  of  the  proceeding  two-dimensional  data  for  use  in  axisymmetric  flows 
past  bodies  of  revolution  may  be  accomplished  by  means  of  Mangler’s  relations^*^*  for  trans- 
forming the  equation  of  n t./on  of  two-dimensional  boundary  layers  to  those  of  equivalent  axi- 
symmetric boundary  layers  on  bodies  of  revolution.  Mangler’s  transformation  relations  are 

f ^ 

djr  . c*  dx 

y - y 

L [42] 

$ - c-^d 
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U(i)  ^ u (») 


Difference  Between  Reynolds  Numbers  of  Neutral  Stobility  and  Tronsition 
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Figure  4 - Position  of  Self-Excited  Transition  as  Function  of  Average  Pressure  Gradient 
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wher«  the  tildes  over  the  quantities  refer  to  equivalent  two*diniensional  conditions  and  C is 

an  arbitrary  constant.  — 

A simple  expression  will  now  be  derived  for  ^ which  can  be  used  to  analyse  two- 

V aSr 

dimensional  data  and  which  can  be  extended  into  an  equivalent  axi  symmetric  system.  First 
the  following  two-dimensional  momentum  equation  is  obtained  from  the  axisymmetric  momentum 
equation,  [80],  by  utilising  Mangler's  relations  [42] 


Since  in  general 

1 d(ifP)  _Cf  dp 

u dSf  V df  V d& 

then  from  [48] 

y ± - JL 

dl  46  6v  df 


[43] 

[44] 

[45] 


Averaging  ^ fiM  over  the  distance  from  neutral  stability  to  transition  in  accordance  with 
V da 

[41]  gives  simply 


i,  di  45  iv 


{0e\, 

*»*  “ */v 


[46] 


In  order  to  use  the  two-dimensional  data  of  Figure  4 to  determine  the  self-excited  tran- 


sition point  for  axisymmetric  flow,  it  is  necessary  to  express  in  terms  of  an  equiva- 

V da 


£.  iU 

V d?' 

lent  axisymmetric  system.  Applying  Mangier’s  transformations  to  [45]  and  integrating  in  ac- 
cordance with  [41]  gives'^  4^  equivalent  axisymmetric  quantities  as 

V da 


J.-JL 

V da  45  5v 


''•i.  - ‘'d 


r 

•/ 

*/v 


da 


[47] 


Furthermore,  in  order  to  use  the  two-dimensional  data  of  Figure  4 to  determine  the  self- 
excited  transition  points  for  axisymmetric  flows,  it  is  necessary  to  convert  the  Reynolds 
number  for  self-excitod  transition  in  two-dimensional  flows  to  that  for  axisymmetric 
flows  on  bodies  of  revolution  Now  Mangler’s  transformations  [42]  give  for  all  Reynolds 

numbers  of  laminar  flows 
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Rq  1481 

Owing  to  the  presence  of  the  arbitrary  constant  C,  it  is  obvious  that  Mangler's  transformations 
are  not  sufficient  to  determine  the  critical  values  of  the  Reynolds  numbers  for  transition. 
Hence,  other  considerations  must  apply.  In  the  analysis  of  the  stability  of  laminar  flows  on 
bodies  of  revolution,  it  was  found*' that  the  resulting  linearised  equation  of  the  disturbed 
flow  at  the  limiting  condition  of  large  Reynolds  numbers  was  the  same  as  that  for  two*< 
dimensional  flows.  Hence 


and  as  stated  previously,  the  neutral  stability  chart  of  Figure  8 is  applicable  to  both  two* 
dimensional  and  axisymmetric  flows  on  bodies  df  revolution.  As  a first  approximation  a 
similar  condition  will  be  assumed  to  hold  for  self-excited  transition.  Then 


[50] 


As  a check  on  the  essential  validity  of  the  preceding,  the  self-excited  transition  points 
were  computed  for  Lyon's^^  bodies  of  revolution  on  the  basis  of  Figure  4.  Owing  to  the  non- 
existence of  low-turbulence  wind  tunnels  at  the  time  of  Miss  Lyon's  tests,  the  condition  of 
least  degree  of  turbulence  represented  in  her  tests  was  that  when  the  turbulence-producing 
screens  were  not  inserted  in  the  wind  tunnel.  Comparison  is  then  made  in  Table  1 of  the  axial 
locations  l/L  of  the  computed  neutral  stability  points  and  self-excited  transition  points  with 
the  test  data  for  the  no-screen  condition. 


TABLE  1 


Axial  Locations  l/L  of  Transition  Data  for  Lyon's  Bodies  of  Revolutions^ 


Subject 

Model  A 

/T^- 2.09x10* 

Model  B 

/?^-  2.075xM* 

Conputed  Neutral  Stability  Point 

0.20 

0.13 

Conpated  Salf-Excitad  Tranaition  Point 

0.S6 

0.7 

Hoaawid  Tranaition  Rtfion 
(without  screen) 

O.SO-0.70 

0.30-0.35 

Regions  of  Favorable  Pressure  Gradient 

0-0.30 

0-0.13 

0.40-0.61 

Reiions  of  Adverse  Pressure  Gradient 

0.30-1 

0.13-0.40 

0.61-1 
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The  measured  transition  points  should  be  behind  the  neutral  stability  points  and 
forward  of  the  self>excited  transition  points  owing  to  the  small  amount  of  natural  turbulence 
in  the  wind  tunnel.  Examination  of  Table  1 confirms  this.  The  reason  for  the  large  difference 
between  the  measured  transition  point  and  the  self*excited  transition  point  on  Model  B is  due 
to  the  favorable  pressure  gradient,  0.40  <;  l/L  < 0.61,  which  tends  to  delay  the  development 
of  selNexcited  transition. 

In  order  to  obtain  the  neutral  stability  and  transition  points,  it  is  more  convenient  to 
use  nondimensional  relations.  From  [36]  there  results 


It  also  follows  that 


[511 


[52] 


153] 


The  neutral  stability  point  or  self*excited  transition  point  is  determined  graphically  by  the 
intersection  of  a curve  representing  the  appropriate  parameters  of  the  given  boundary  layer 
with  the  curves  of  Figure  3 or  4,  as  the  case  may  bo. 

The  position  of  transition  on  a body  in  a wind  tunnel  or  other  flow  facility  depends, 
to  a great  extent,  upon  the  scale  and  intensity  of  the  turbulence  in  the  free  stream.  To 
simplify  the  analysis  the  representative  measure  of  turbulence  in  the  stream  will  be  taken  as 
the  ratio  of  the  root  mean  square  of  the  velocity  fluctuations  in  the  ^-direction  to  the  mean 
velocity  of  the  f^  stream  or  Data  for  flat  plates  without  pressure  gradient  from 


tests  in  wind  tunnels  are  plotted  in  Figure  5,  the  ordinates  being  the  difference  in 

the  boundary-layer  Reynolds  numbers  for  transition  and  that  for  neutral  stability 

and  the  abscissae  being  the  percentage  of  turbulence  in  the  free  stream.  It  is  to 

be  noted  that  for  flat  plates  the  neutral  stability  point  seems  to  be  the  limiting  position  of 
transition  for  streams  of  increasing  turbulence  while  the  self-excited  transition  point  is  the 
limiting  position  of  transition  for  streams  of  decreasing  turbulence. 
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A similar  trend  for  boundary  layers  with  pressure  gradients  is  indicated  by  the  test  data 
for  the  flow  over  an  elliptic  cylinder*^  as  plotted  in  Figure  5.  On  the  basis  of  the  neutral 
stability  point  and  the  self-excited  transition  point  as  the  lower  and  upper  limiting  positions 
of  transition  for  boundary  layers  with  pressure  gradients  in  general,  a rough  criterion  may  be 
established  for  estimating  the  position  of  transition  on  bodies  immersed  in  streams  containing 


turbulence.  If  the  curve  for  transition  for  any  boundary  layer  with  a pressure  gradient  is  as* 
sumed  geometrically  similar  to  that  for  a flat  plate  as  plotted  in  Figure  5,  the  Reynolds 
number  for  transition  for  boundary  layers  on  bodies  in  a stream  with  turbulence  may  be 
stated  as 

^6,tb  ~ , / ^0,  tb  - ^d,N\  [54] 

^e.ts  ~ ^d.N  ~ ^d.Nf 

0 

where  the  subscript  0 refers  to  flat*plate  values  without  pressure  gradient.  For  more  precise 
results,  it  is  best  to  depend  upon  measured  locations  of  the  position  of  transition  to  obtain 
values  of  the  Reynolds  numbers  for  transition  characteristic  of  the  particular  wind  tunnel  or 
flow  facility  being  utilized. 

AXISYMMETRIC  TURBULENT  BOUNDARY  LAYER 


GENERAL 

Owing  to  the  incomplete  state  of  present  knowledge,  concerning  the  mechanics  of 
turbulent  flow  processes,  the  analysis  of  turbulent  boundary  layers  lacks  the  clearly  defined 
features  of  laminar  boundary  layers.  In  order  to  arrive  at  results  of  immediate  utility,  oxten* 
sive  reliance  has  to  be  placed  on  empirical  data  to  augment  theoretically  derived  relations. 

Turbulent  flows,  in  general,  may  be  treated  from  the  Reynolds  viewpoint,  this  consid* 
ers  turbulent  flow  to  consist  of  a mean  flow  upon  which  a fluctuation  flow  of  much  smaller 
magnitude  is  superimposed.  After  the  combined  mean  and  fluctuation  quantities  are  substi- 
tuted into  the  Navier-Stokes  equations  of  motion  for  viscous  flow,  appropriate  time  averages 
of  the  resulting  flow  lead  to  the  Reynolds  equations  of  motion  containing  separate  terms  for 
the  mean  quantities  and  the  fluctuation  quantities.  The  form  of  the  Reynolds  equations  is 
similar  to  the  Navier-Stokes  equations  with  the  significant  exception  of  the  presence  of  addi- 
tional terms  which  consist  of  averages  of  various  products  of  the  fluctuation  velocities. 
These  additional  terms  may  be  shown  to  act  as  apparent  stresses  (Reynolds  stresses)  within 
the  flow.  It  is,  however,  ihe  present  lack  of  analytical  relations  for  the  Reynolds  stresses 
which  has  made  the  theoretical  treatment  of  turbulent  flows  so  difficult. 

Turbulent  boundary-layer  equations  are  derived  from  the  Reynolds  equations  in  the 
same  way  as  laminar  boundary-layer  equations  are  derived  from  the  Navier-Stokes  equations, 
i.e.,  by  the  process  of  eliminating  terms  of  negligible  magnitude. 

In  the  case  of  axisymmetric  flow  past  a body  of  revolution  with  negligible  longitudinal 
curvature,  the  Reynolds  equations  of  motion  become* 


*Tb»  Rayaolda  •quatlona  of  motion  for  aonorol  currUlnear  coordloatoo  ara  wrlttaa  In  tanaor  notation  in  Refar> 
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u + l -if  JL  Fr/iif- + J.  ^ 

d»  dy  p d»  rdy  Uy  d»/J  pf  5*  L V /J 

. j.  ± r,/-^:v)],2! 

/J  f dx 

u ±s.  + ui»  ^ 1 i£  .±  ± fff^  - ~ — F'  f-  P^)\ 

dy  By  p dy  r Bx  ^ 5y/J  p'  ^y  L V /J 
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and  the  equations  of  continuity 


iiLH2  . iiHil  - 0 

Bx  ay 


[56] 

[57] 


The  unprimed  quantities  refer  to  the  mean  flow  and  the  primed  quantities  to  the  fluctuation 
flow.  The  bars  over  the  various  products  indicate  averages  in  accordance  with  the  Reynolds 
concept.  Specifically,  u,  v,  primed  or  unprimed,  represent  components  of  velocity  in  the 
X,  y,  and  <f>  directions,  respectively.  As  before,  x is  the  distance  along  the  meridian  profile 
and  y is  the  distance  perpendicular  to  the  profile,  tft  is  the  angle  between  any  meridian  plane 
and  some  reference  meridian  plane. 

Eliminating  terms  of  negligible  magnitude  in  [55]  gives  the  following  equations  of 
motion  for  the  axisymmetric  boundary  layer  on  a body  of  revolution 

^O'P)  + J_  ^(>‘^) 

Bx  By  p Bx  or  Bx  or  By 


o 

■ 

[58] 

9 

where  —r 

a - pv'* 

is  the  negative  value  of  a Reynolds  normal  stress  and 

[59] 

T - p ij£  - pu'v" 

[60] 

is  the  total  shearing  stress. 

Integrating  the  boundary 'layer  equations,  [58],  across  the  boundary  layer  in  the  y- 
direction  from  y - 0 to  y - 5 and  incorporating  the  equation  of  continuity,  [56],  gives  the  dif> 
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ferential  form  of  the  momentum  equation  for  axisymmetric  turbulent  boundary  layers  on  bodies 
of  revolution 


+ (A  ♦ 2) 

d9 


i ilL 

V dn 


[61] 


where  the  momentum  area  Q is  deHned  in  Equation  [82].  In  terms  of  the  displacement  area  A , 
defined  in  [10]  and  fl,  the  shape  parameter  for  axisymmetric  flow  is 


A - -A-  [62] 

Q 

The  momentum  equation,  [61],  will  now  be  considered  for  the  two  cases  where  the  boundary 
layer  is  thin  compared  to  the  radius  of  the  body  and  where  the  boundary  layer  is  of  like 
magnitude. 


BOUNDARY  LAYER  THIN  RELATIVE  TO  BODY  RADIUS:  £ « 

In  the  region  forward  of  the  tail,  the  boundary  layer  is  thin  relative  to  the  radius  of  the 
body  B «r^.  Hero  the  momentum  equation,  [61],  may  bo  reduced  to  a simpler  form. 

Since  within  the  boundary  layer 


[63] 

it  follows  from  [32]  that 

^6  < 0 

< K * ^ 

[34] 

from  [10]  that 

r^B*<  A* 

< + 5) 

[64] 

and  from  [62] 
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where  H ie  the  two-dimensional  shape  parameter 
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Evidently  for  8 « , fi  - 
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Furthermore 


S 8 

f I s j »■'» 


[68] 


since 


/ f * \ dr 

adyj-^  is  empirically  of  a lower  order  of  magnitude.  Hence  the  momentum 


equation,  [61],  becomes  for  S « r^ 

. (H  . 2)  Isif  . r fji  » -1-  ^ f orf,] 

' ’ (/  d,  J n 
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[69] 


T C ^ 

of  terms  — and  ^ x — | ady  may  be  considered  as  an  effective 
pV*  p dx  \ 

...  Jo 


The  sum 

coefFicient  of  local  skin  friction,  i.e.. 
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o dy 
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In  their  study  of  two-dimensional  boundary  layers  in  adverse  pressure  gradients,  Wieghardt 
and  Tillmann  found  that  the  values  of  ( ''  P^^)tu  computed  from  the  momentum  changes  of 

the  mean  flow  displayed  unexpected  increases  near  separation.^*  Since  tests  by  Ludwieg  and 
Tillmann^^  have  shown  the  values  of  / pU‘  to  decrease  in  an  adverse  pressure  gradient 
faster  than  those  for  flat  plates  without  pressure  gradients,  the  explanation  of  the  increase  in 
( T,^/p  has  to  be  sought  elsewhere.  Wieghardt  has  attributed  this  apparent  increase 

wholly  to  the  convergence  of  the  flow  caused  by  the  thickening  of  the  boundary  layer  on  the 
opposite  sides  of  the  wind  tunnel. On  the  other  hand,  Newman^^  and  various  investigators 
of  the  National  Advisory  Committee  for  Aeronautics^® have  shown  this  increase  to  be 

1 J /•  5 

partly  accounted  for  by  the  increase  in  the  value  of  the  normal-stress  term  -=—x  — I ody, 

pV  dx  Jq 

especially  close  to  separation.  Hence  for  adverse  pressure  gradients,  the  decrease  in  the 
value  of  r..  / pU^  tends  to  be  compensated  by  the  increase  in  the  value  of 


d 

Hx 


I 


w p 
S 


ady.  This  makes  the  coefficient  of  local  skin  friction  for  flat  plates  without 


JO 

pressure  gradient  ( p^^)q  ^ close  approximation  to  the  effective  skin  friction,  at  least 

for  moderate  pressure  gradients.  Lyon’s  experimental  results  tend  to  substantiate  this.^^ 
Therefore  with 


(M  ■ (>) 


[71] 


«ir 


Equation  [69]  becomes 


Coefficient  of  Locol  Skin  Friction  for  Flat  Platt* 
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+ («  + 2) 
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ilL 

dx 
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The  coefficient  of  local  skin  fricUon  for  flat  plates  (.  r^/  pLI^)^  decreases  slowly 
with  local  Reynolds  number  Rq.  A suitable  formula  valid  for  a large  range  of  Reynolds  num- 
bers is  that  derived  in  Reference  38  from  the  Kannan-Schoenherr  formula  for  the  total  dr^j;  of 
flat  plates, 


0.01466 


2 ■♦•0.4348 


which  is  plotted  in  Figure  6. 


[73] 


A rapid  procedure  (and  one  sufficiently  accurate  for  drag  calculations)  is  to 
approximate  {r^/  pV^)^  with  a power  law  and  to  assume  H constant.  With 


<0 


[74] 
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where  m and  are  constants,  [72]  becomes 


<h 


(//  + 2)  dU  , 
* — 


U 


[75] 


which  is  a flrst^order  differential  equation  in  r^  $ of  the  Bernoulli  type.  For  H constant,  [75] 
gives  for  interval  of  integration  to 


where 


1 + 1 
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i-m 


dx 


[76] 


}.  (!  + «)(« + 2)  [77] 

Lyon  found  that  H - 1.4  gave  good  agreement  between  the  experimental  and  calculated 
drags. Since  (r^/  pU^)^  in  Equation  [78]  gives  almost  a straight  line  on  a log-log  plot 
as  shown  in  Figure  6,  a power  law  approximation  proves  a close  fit.  For  10^  < R0<  10*,  a 
least-squares  fit  of  [73]  gives 


at  s 0 • 16B6 

-0.008361 


[78] 


A more  convenient  nondimensional  form  of  [76]  is 


i^ysir):  -{mi 


/I  V ^ 

(1  + ■•)  Co 
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BOUNDARY  LAYER  NEAR  THE  TAIL 


It  is  apparent  that  the  previous  assumption  5 « r^  is  no  longer  valid  towards  the  after 
end  of  the  body  where  the  radius  goes  to  sero  at  the  tail  and  the  boundary  layer  progres- 
sively thickens.  Consequently  the  general  form  of  the  momentum  equation  [61]  has  to  be  em- 
ployed; for  constant  A,  it  may  be  integrated  to  the  form 


+ ±-  ± 

pU^  ^ 


dx 

[80] 


where  the  subscript  g refers  to  the  initial  point  of  integration  and  the  subscript  e to  the  tail 

of  the  body.  ^ 

Experimental  evidence indicates  that  r -lllL  •».  -L_  A.  f atdy  in  the  rights 

pl/2  pV^  dx 

hand  side  of  Equation  [80]  is  substantially  linear  with  respect  to  x along  the  after  end,  and 
drops  to  Eero  at  the  tail.  Consequently 
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F nrth<«t7nom,  lh(«  vArittion  of  V with  raspwct  to  « is  also  substantially  linear  along  the  after 

onH  or 


U - - X 


[82] 


Inserting  theae  linear  relations,  Equations  [81]  and  [82],  into  momentum  equation  [80]  gives 


^(A  Vy*  - (A  + 4) 


[83] 


or  aoadimensionally 


t U»  iHitlal  |H«Milion  y of  Uie  tail  portion  of  the  boundary  layer  may  be  taken  at  0.8  of  the 
iviuiik  Iht*  iHMiy.  An  average  value  of  A « 1.4  over  the  tail  portion  is  appropriate  in  view  of 

aUo  nf  A - 1.42  found  on  the  Mark  13  Torpedo.^ 

In  (he  region  of  the  boundary  layer  near  the  tail,  it  is  useful  to  establish  relations  be* 
« ami  (>,  3*  and  A*,  and  A and  H.  A power  law  is  assumed  for  the  velocity  profile  or 


u 

U 


[85] 


)ub.)Ututiug  thiH  expression  into  Equation  [32]  and  performing  the  indicated  integration 


lu  III 
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Also  from  Equation  [10] 


[861 


Furthermore, 
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In  the  wake  where  ~ 0 
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TURBULENT  WAKE 

The  boundary  layer  leaving  the  tail  of  the  body  as  the  wake  has  a somewhat  higher 
pressure  and  larger  momentum  area  0,  than  the  wake  far  downstream.  Since  the  momentum 
area  of  ihe  wake  far  downstream  determines  the  drag  of  the  body,  a relation  is  required 
between  and  the  momentum  area  fl^  of  the  boundary  layer  at  the  tail. 

With  no  skin  friction  in  the  wake  and  with  negligible  effect  from  the  normal-stress  term 
the  momentum  equation  [61]  reduces  to 


^ + (A  + 2)il^-0 

d*  ' V dx 


[90] 


Integrating  by  parts  over  the  length  of  the  wake  from  the  tail  to  infinity  downstream  yields 

k. 


r r * 


[91] 


where  the  limiting  value  of  A at  infinity  is  unity. 
The 


r " 1 

evaluation  of  exp  I In  — dA  proceeds  empirically.  Experimental 

Ui  ^ J 


evidence 


suggests  an  empirical  fit  of  a higher  order  parabola  of  form 


In 


JL.  . 

U.-i/ 


u. 
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Figure  7 • Variation  of  Axisymmetric  Shape  Parameter  h in  Wake 


Then 


and 


Finally 


(ft,  + 2)  9 + 3 


[93] 


[94] 


[95] 


This  completes  the  relationships  required  for  calculating  the  drag  from  the  boundary-layer 
development.  It  is  seen  that  both  the  momentum  area  and  the  shape  parameter  at  the  tail  are 
required  for  a precise  solution  of  the  drag. 

The  choice  of  a value  for  q is  not  critical  since  a large  difference  in  q results  in  only 
a small  error  in  the  drag.  Considering  only  a very  limited  amount  of  two-dimensional  data, 
Young^  uses  a linear  relation  for  Equation  [92],  that  is,  a value  of  q of  unity.  As  shown  in 
Figure  7,  test  data  for  a body  of  revolution,  the  Mark  18  Torpedo,^^  though  meager,  indicate 


that  the  ordinate 


In  V^/V 


approaches  the  limit  of  zero  at  the  origin  more  gradually  than  by 


a linear  relationship.  Until  additional  test  data  are  available,  it  is  suggested  that  a value  of 
q ml  he  used  for  drag  calculations.  If  the  value  of  - 1.42  (from  the  data  for  the  Mark  13 
torpedo)  is  also  used  in  Equation  [95],  there  results  finally 


n. 


[96] 


CONCLUDING  REMARKS 

The  method  described  in  this  report  gives  a relatively  simple  procedure  for  computing 
the  viscous  drag  of  a body  of  revolution  to  an  engineering  degree  of  accuracy.  Any  marked 
deviation  between  the  calculated  and  measured  values  of  drag  should  be  expected  only  in  the 
case  of  a shape  differing  radically  from  the  usual  streamlined  figure. 

To  test  the  method  of  this  report,  drag  coefficients  of  the  model  of  tlie  airship  AKRON 
were  computed  at  various  Reynolds  numbers  for  comparison  with  measured  values. 

Measured  values  of  the  pressure  distribution  and  measured  locations  of  the  transition  point 
were  used  in  the  calculations.  The  computed  values  were  found  to  lie  between  the  measured 
values  as  shown  in  the  accompanying  table. 

TABLE  2 


Comparison  of  Measured  and  Computed  Values  of  Drag  Coefficient  of 
1/40-Scale  Model  of  Airship  AKRON 


Reynolds 

Number 

Drag  Coefficient  Cp 

Experimental 

Computed 

Wooden 

Model 

l2.3x  # 

0.0198 

0.0228 

15.0  X 10* 

0.0193 

0.0223 

17.  3 X 10* 

0.0190 

0.0219 

0.0211 

In  conclusion,  it  is  noted  that  in  the  case  of  low  Reynolds  numbers  the  viscous  drag  of 
bodies  of  revolution  depends  to  a large  measure  on  the  position  of  transition,  and  that  in  the 
case  of  high  Reynolds  numbers,  it  depends  on  the  detailed  development  of  the  turbulent  bound' 
ary  layer  in  pressure  gradients.  The  current  theories  on  turbulent  flow  are  semi-empirical  in 
nature.  Accordingly,  Before  all  the  processes  determining  viscous  drag  are  fully  elucidated, 
there  is  need  for  accurate  measurements  of  such  factors  as  velocity  profiles  and  shearing- 
stress  profiles  in  turbulent  boundary  layers  especially  at  high  Reynolds  numbers. 
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For  convenient  reference,  the  various  steps  involved  in  calculating  the  viscous  drag  are 
summarized  in  Appendix  I. 
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APPENDIX  1 

CALCULATION  PROCEDURE 

For  convenient  reference,  the  various  steps  for  computing  the  viscous  drag  of  bodies 
of  revolution  in  axisymmetric  flow  are  listed  here. 

1.  Prescribed  Data 


a.  Profile  dimensions 


b.  Pressure  distribution 


P-Pm 
2^  - 


■(tr  •'«) 


c.  Body  Reynolds  number 


2.  Transition  Point 


N..U.UUW.U,  point  (i)^ 

(1)  He-  i)  is  computed  from 


[511 


is  computed  from 


► 4) 


<*> ■ '{t  S)  is  plotted  to  intersect  the  neutral  stability  curve  of  Figure  3 

U,giv,(i)^ 

b.  Self-excited  transition  point 


88 


(1)  ^ fli\  is  computed  for  various  assumed  values  of/i- ) from 

V dV  \Llf,  \Ll,^ 


? iS’  - i_  li 

V d*  45  5 ' 


{!iWL 


(t:) 


U L 


UJ 


W) 


N 


I 


[53] 


(2) 


Rq  Rq  ^.jis  plotted  to  intersect  the  curve  of  Figure  4 to  give 


correct  value  of  LLi 
Wts 


Ce  Transition  point  for  free-stream  iurbulencej 

(1)  For  prescribed -!L — of  free  sUeam  is  determined  from 

Figure  5. 

(2)  fie,, 6 determined  from 

fie.tb  " f^$,N  _ / ~ ^,jv\ 

fie.u  " ^e.N  “ f^e.Nfo 

3.  Laminar  Boundary  Layer 

at  transition  point^-^^  is  computed  from 


[54] 


(^)!  i’Jf.  ■ « 1, 

' " *«  ■ '(f) 


/r  \2 


(t)' 


has  been  computed 


U'J,  \ul  Wl.lHA 
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4.  Turbulent  Boundary  Layer 

a.  Thin  boundary  layer,  8 «r^ 

Q is  computed  over  distance |-i>]  0* 


8 from 


(s):  ■amt 


1 \/-»  jfl\ 

[»• 

(</Mi 

t) 

J secad(j-J 

where  m - 0.1686 
; - 3.9732 
- 0.006361 

b.  Boundary  layer  near  tail 
(1)  ( is  computed  from 

\L^L 


(a  fe). 


ra.'ftKSr 

- ii\  sec  a 

±)  5 

(“•'<** Hit)  (it -it) 


A value  of  A ■ 1.4  and^j 


0.8  is  to  be  used. 


is  computed  from 


6.  Viscous  Dreg 

Dng  coefficient  Cq  is  computed  from 
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